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. P. ; Abstract 

A theoretical model of high-Tc Josephson Field Effect Transistors (JoFETs) 
based on a Ginzburg-Landau free energy expression whose parameters are 
field- and spatially- dependent is developed. This model is used to explain ex- 
perimental data on JoFETs made by the hole - overdoped Ca-SBCO bicrystal 
junctions (three terminal devices). The measurements showed a large modu- 
lation of the critical current as a function of the applied voltage due to charge 
modulation in the bicrystal junction. The experimental data agree with the 
solutions of the theoretical model. This provides an explanation of the large 



field effect, based on the strong suppression of the carrier density near the 
grain boundary junction in the absence of applied field, and the subsequent 
modulation of the density by the field. 
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The development of practical high-Tc three terminal devices has received much attention 
recently. There are many different directions that have been investigated: the electric field 
effect transistor, the fiux-fiow transistor, quasiparticle injection devices, etc. The super- 
conducting field effect transistors (suFETs) with a homogeneous thin film channel [0,0 and 
JoFETs with a Josephson junction channel are electric field related devices. In general, 
the required electric field for the field effect in JoFETs is orders of magnitude smaller than 
the field required in suFETs with homogeneous thin film. 

A basic question that therefore arises is to understand the large effect that an electric 
field has on the transport properties in JoFETs. Dong et al. showed that there is a 
23% modulation of the critical current on a grain boundary of a 50 nm thick channel of 
Smi_a;Caa;Ba2 Cn^Oy (Ca-SBCO). Nakajima et al. reported a 5% modulation of critical 
current on a 60 nm thick YBCO grain boundary junction channel which can be described 
by a parallel resistor model. A modulation of several per cent (maximum 8% in Ic) has 
been reported by Petersen et al. ||^ on the transport properties of a less than 32 nm thick 
YBCO grain boundary junction channel both in the normal and superconducting state. 
The most recent JoFET experiments were carried out by Mannhart's group where a 
similar field modulation of Ic was observed. However it is still not clear why all suFETs 
with grain boundary junction (JoFETs) show much bigger effects than suFETs made of 
homogeneous film under the same applied field. Candidate explanations include a weakly 
coupled SNS model in the dirty limit 0], a parallel resistor model at high bias current, and 
the electromechanical effect in the dielectric layer at low bias currents 0. But none of the 
above gives a quantitative explanation of the dependence of the critical current on the field. 
The mechanism responsible for this large field effect therefore remains an open and very 
important question. 

In this paper we attempt to clarify this issue, using a phenomenological model based on 
the Ginzburg-Landau (GL) theory of phase transitions. A similar model has been already 
developed for the case of YBa2Cu307_5 grain boundaries in bicrystals 0. The basic result 
of that study was that the oxygen depletion can account for a major portion of the change 
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from weak to strong coupling of grain boundaries, which is experimentally observed |T0 
as the misorientation angle is increased. The modification of the oxygen content leads to 
the variation of the critical temperature as a function of distance from the boundary. The 
detailed way that this occurs is not well understood from a microscopic point of view. The 
phenomenological approach of Ref. allows one to simulate the behavior of the system in 
terms of a few measured parameters and to calculate electromagnetic properties. Our aim 
in this paper is to apply this method to the field efffect modulations seen in the JoFETs 
of Ref. 13,1] . We demonstrate that the JoFET systems can be fit into the same conceptual 
framework as the YBa2Cu307_5 grain boundaries. 

The experimental measurements 0] were taken on hole-overdoped Ca-SBCO bicrystal 
junctions with 30% doping of Ca (x = 0.3) at 20 K and 4.2 K. The junction itself is a 
grain boundary with a misorientation of 24°, created by growing the high-T^ film on an 
SrTiOs substrate with such a boundary. A schematic view of the experimental setup is 
shown in Fig. 1. Junctions of this kind in the YBa2Cu307_5 system have been shown to be 
oxygen-deficient []Tn|,|n| with a consequent lowering of the critical temperature. 

The basic experimental result is the modulation of the (normalized) critical current as a 
function of the applied gate voltage, as shown in Fig. 2. The input needed for the theory, 
is the critical temperature Tc as a function of the distance from the boundary plane, and as 
a function of the applied field. The theory then gives a prediction for the critical current. 



Detailed work on the effect of doping on the critical temperature of the cuprates |1T2|JT3 
shows that Tc follows a parabolic relation: T/Tc^^ax = 1 ~ 82.6 * {n — 0.16)^, where n is the 
carrier concentration (holes per Cu02 unit) and T^^max = 70K for the material under study 
(Cao.3Smo.7Ba2Cu30y) 0. From the specific characteristics of the specimen, given that a 
grain boundary of 24° corresponds to a weakly-coupled bicrystal, then the dependence of 
the concentration n on the distance follows an exponential function in accordance with the 
previous work So at zero applied gate voltage : n{x) = 0.210 — 0.206 * exp(— 0.22;/^), 
where ^ is the superconducting coherence length which has an approximate value of 2nm. 
Also, the effect of the applied gate voltage Vg is taken as a linear contribution in the 



concentration function, since the induced charge density AN for the specific material can 
be found from the observation that AN/Vg = Cg/\e\ ~ 3.2 * 10^^cm~'^V~^ where Cg is 
the areal gate capacitance and e is the electron charge. The contribution is then 6n = 

In the one-dimensional case we consider, the GL free energy that has to be minimized 
takes the form (we use a gauge where the vector potential A = for convenience, since no 
magnetic field is present and the free energy is gauge invariant): 

F = a{x,T)^' + ^^'-^{^^r, (1) 
2 2m* ax 

where m* is the effective mass. We take the following form for a{x,T): 

a{x,T) = aoTe(x)'[tanh(3.0v/T,(x)/T - 1)]^ ifT,>T (2) 
= QaoT\T,{x)/T - 1) if T, < T. (3) 

This expression for a{x, T) for T < Tc is an analytic fit to the strong-coupling form of the 
gap function in Bardeen-Cooper-Schrieffer theory. Above we have less knowledge about 
the form of the coefficients of of GL theory. We have chosen a form for a in this regime 
guided by two considerations: the expression for a should continue smoothly through Tc, 
and should be consistent with the equation ^ = (^^/2m*a)^/^, where ^ is the coherence 
length. A function different from the usual first order term (T — Tc) in the expansion close 
to Tc is necessary in order to cover the whole range of temperatures (it is stressed here that 
Tc varies with the distance). It is obvious that if we expand the above function close to 
Tc we get the usual (T — Tc) term. (3 is taken as constant. The theory thus involves the 
assumption that GL theory may be used over a broad range of temperatures. The great 
virtue of using this approach instead of a microscopic theory is that the parameters can 
be related to a number of observable quantities. The above choice has been tested in the 
succesful fitting or prediction of several quantities (calculation of the order parameter, NMR 
studies, specific heat etc.) [jl^,|T5[. 

It is now convenient to write: \I' = |\l/|exp(i0) in which case the current density J is 
given by: J = (he/m*) * |\E'p * d(f)/dx 
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These expressions may be simphfied by the definitions: f{x) = \I//\I/(oo) , h{x) = 
a{x,T)/a{oo,T) and j = J/ Jc{oo), while x is taken in units of ^. 
Here Jc(oo) is the bulk depairing current : 

J,(oo) = 2e|vl^(oo)P^(^^^^)V^. (4) 

The Euler-Lagrange equation corresponding to the GL free energy in Eq. |l] is the differ- 
ential equation: 

^ - (4jV27f ) + h{x)f -f = 0, (5) 

The computational problem is to solve the nonlinear differential equation (13) with the 
boundary conditions: /(±oo) = y/y and /'(±oo) = , where y is the solution of the equation 
■.y'-y' = (4/27)j2. 

At low current densities j, a superconducting solution exists. The critical current is 
found by increasing j until no superconducting solution exists. This yields the quantity 
Ale/ Ico where Ale = hiVg) — -^c(O) = IciVg) — ho- We plot the calculated results together 
with the experimental data in Fig. 2. It's clear that the model under study can reproduce 
the basic experimental finding of the large electric field effect on grain boundary Josephson 
junctions. The results show modulation in critical current from a few percent for small 
voltages to almost 20% for the largest values that have been used. 

The calculation of the order parameter gives a very good picture of the suppression 
of superconductivity as we approach the boundary (Fig. 3). Interestingly, there is also a 
region with order parameter values l^'l greater than the bulk value of \E' due to the shape 
of the function Tc{x), leading to an enhancement of the superconductivity just before the 
surpression. Temperature plays a more important role in the shape of the order parameter, so 
at higher experimental temperature there is a well defined nonsuperconducting region close 
to the grain boundary. The application of the gate voltage alters the critical current which 
has small effect on the order parameter (due to the fact that it is several order of magnitude 
less than the depairing current at infinity). This small effect can be barely detected in 
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the asymptotic values of \1/ as well as the value of the peak of Furthermore, similar 
calculations in the underdoped regime (e.g. in an as-made Y Ba2Cu^0y film on bicrystal) 
have been performed and the quantity AIc/Ico shows that this case is less sensitive to the 
field (about half the modification observed in the overdoped regime). 

The calculations demonstrate that if there is relatively weak coupling between the two 
sides of the boundary, a relatively small modification of the carrier density can have dramatic 
consequences. The boundary serves effectively as a proximity-effect junction which changes 
from S — S' — S towards S — N — S as the field is applied. The result is the observed large 
field effect. Thin films do not show the same effect because the field is being applied to a 
strong superconducting region. 

Another point to emphasize is that the above model doesn't distinguish between s-wave 
or d-wave symmetry of the gap function and consequently of the order parameter. The GL 
theory has the identical form for the two cases [|l^]. Thus the actual microscopic mechanism 
of high-Tc superconductivity doesn't affect our results. 

These results provide guidance for further investigations in this field. The large field 
effect in JoFETs can be accurately predicted within the Ginzburg-Landau theory while the 
complications of the microscopic theory are avoided. 

This work was supported by the NSF Grant. No. DMR-9214707 (J.B. and R.J.), by 
NSF Grant. No. DMR-9404579 (Z.-W.D. and T. V.) and by the Dutch Program for High 
Temperature Superconductivity (P.H.). 
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FIGURES 

FIG. 1. Schematic view of JoFET and the circuit outhne for the field effect measurements 

FIG. 2. In Fig.2(a) the experimental data from Ref.[2] (the circles represent the 15 /xm wide 
device and the squares the 30 /xm wide device), in Fig.2(b) the results of the calculations from 
(GL) theory for the two experimental temperatures T=20K (solid line) and T=4.2 K (dashed line) 

FIG. 3. The calculated order parameter for the gate voltage of 40 Volts and for the two 
experimental temperatures 4.2 K and 20 K 
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